The diurnal nature of future extreme precipitation intensification

Short‐duration, high‐impact precipitation events in the extratropics are invariably convective in nature, typically occur during the summer, and are projected to intensify under climate change. The occurrence of convective precipitation is strongly regulated by the diurnal convective cycle, peaking in the late afternoon. Here we perform very high resolution (convection‐permitting) regional climate model simulations to study the scaling of extreme precipitation under climate change across the diurnal cycle. We show that the future intensification of extreme precipitation has a strong diurnal signal and that intraday scaling far in excess of overall scaling, and indeed thermodynamic expectations, is possible. We additionally show that, under a strong climate change scenario, the probability maximum for the occurrence of heavy to extreme precipitation may shift from late afternoon to the overnight/morning period. We further identify the thermodynamic and dynamic mechanisms which modify future extreme environments, explaining both the future scaling's diurnal signal and departure from thermodynamic expectations

On skillful decadal predictions of the subpolar North Atlantic

The North Atlantic is a crucial region for the prediction of weather and climate of North America and Europe and is the focus of this analysis. A skillful decadal prediction of the surface temperature was shown for several Earth system models, with the North Atlantic standing out as one region with higher predictive skill. This skill assessment concentrates on the rapid increase of the annual mean sea surface temperature of the North Atlantic subpolar gyre by about 1 K in the mid‑1990s and the adjacent years. This event-oriented analysis adds creditability to the decadal predictions and reveals the potential for improvements. The ability to simulate the observed sea surface temperature in the North Atlantic is quantified by using four versions of decadal predictions, which differ in model resolution, initialization technique, and the reanalysis data used in the assimilation run. While all four versions can reproduce the mid-1990s warming of the subpolar North Atlantic, the characteristics differ with lead time and version. The higher vertical resolution in the atmosphere and the higher horizontal resolution in the ocean improve the decadal prediction for longer lead times, and the anomaly initialization outperforms the full-field initialization for short lead times. The effect from the two different ocean reanalysis products on the predictive skill is strongest in the first two prediction years; a substantial cooling instead of the warming in the central North Atlantic reduces the skill score for the North Atlantic sea surface temperature in one version, whereas a too large interannual variability, compared with observations, lowers the skill score in the other version. The cooling patches are critical since the resulting gradients in sea surface temperature and their effect on atmospheric dynamics deviate from observations, and, moreover, hinder the skillful prediction of atmospheric variables

Large-scale secondary circulations in a limited area model – the impact of lateral boundaries and resolution

Within their domain, regional climate and weather forecasting models deviate from the driving data. Small-scale deviations are a desired effect of adding regional details. There are, however, also deviations of the large-scale circulation, which can be caused by orographic effects and depend on the large-scale flow condition. These ‘secondary circulations’ (SCs) are confined to the model domain due to the prescribed boundary conditions. Here, the impact of different regional model configurations on the SC is analysed in a case study for the European region using an ensemble approach. It is shown that at 500 hPa, vortices of the SC have diameters on the order of several thousand kilometres and are related to wind speed anomalies of more than 5 m/s and geopotential height anomalies of more than 5 dam. The spatial structure and the amplitude of the SC strongly depend on the location of the lateral boundaries. The impact of the boundary location on the anomalies is on the same order of magnitude as the anomalies themselves. The resolution of the regional model, as well as the application of spectral nudging and a smoothed topography, affects mainly the amplitude of the SC, but not the spatial structure

Kinematic vorticity number - a tool for estimating vortex sizes and circulations

The influence of extratropical vortices on a global scale is mainly characterised by their size and by the magnitude of their circulation. However, the determination of these properties is still a great challenge since a vortex has no clear delimitations but is part of the flow field itself. In this work, we introduce a kinematic vortex size determination method based on the kinematic vorticity number Wk to atmospheric flows. Wk relates the local rate-of-rotation to the local rate-of-deformation at every point in the field and a vortex core is identified as a simply connected region where the rotation prevails over the deformation. Additionally, considering the sign of vorticity in the extended Wk-method allows to identify highs and lows in different vertical layers of the atmosphere and to study vertical as well as horizontal vortex interactions. We will test the Wk-method in different idealised 2-D (superposition of two lows/low and jet) and real 3-D flow situations (winter storm affecting Europe) and compare the results with traditional methods based on the pressure and the vorticity fields. In comparison to these traditional methods, the Wk-method is able to extract vortex core sizes even in shear-dominated regions that occur frequently in the upper troposphere. Furthermore, statistics of the size and circulation distributions of cyclones will be given. Since the Wk-method identifies vortex cores, the identified radii are subsynoptic with a broad peak around 300-500km at the 1000 hPa level. However, the total circulating area is not only restricted to the core. In general, circulations are in the order of 107m2/s with only a few cyclones in the order of 108m2/s

Application of an object-based verification method to ensemble forecasts of 10‑m wind gusts during winter storms

The object-based method SAL (Structure, Amplitude and Location) was adapted for investigating the errors of forecasts of extreme 10‑m wind gusts associated with winter storms in Germany. It has been applied to a statistically downscaled version of the 51 member ECMWF (European Centre for Medium Range Weather Forecasts) operational ensemble forecast. The horizontal resolution of both downscaled data and of the German weather service's operational analysis data used for verification is 7 km. Forecast errors are subdivided in terms of storm intensity, location and extent. After identifying a set of storm events, objects of moderate and intense 10‑m wind gusts were identified with a local percentile-based threshold (90th percentile for moderate and 98th percentile for intense gust objects). Depending on the intensity of the storm, the gust objects differ in terms of size, shape and intensity. The characteristics of the ensemble forecasts of 10‑m wind gusts can basically be assessed in two different ways. Individual forecast members can be evaluated with respect to the location, intensity and extent of the gust field, and then address the ensemble characteristics by the score distributions. Alternatively, the gust fields' location, intensity and extent can be evaluated by directly using the ensemble mean forecast instead of the individual members. The results of the identified set of storms clearly indicate a high case-to-case variability in the predictability of 10‑m wind gusts objects, particularly when focusing on the structure of intense wind gust objects. It is found, that the gust fields' location and overall intensity can be better estimated from the ensemble mean forecast, compared to the individual forecast members. From a forecaster's perspective this means, that a storms' location and intensity can be well estimated by considering the ensemble mean wind forecasts. Considering the structure of the gust objects, results are different. While for longer lead times, there also seems to be a benefit from applying ensemble averaging, at short lead times the ensemble mean forecast performs equally or worse than most of the individual forecast members. The amplitude error is often the smallest component of the three error types. The findings are particularly relevant when deriving warning information, by giving guidance to forecasters when interpreting ensemble forecasts for severe storms

Weather impacts on various types of road crashes: a quantitative analysis using generalized additive models

Adverse weather conditions can have different effects on different types of road crashes. We quantify the combined effects of traffic volume and meteorological parameters on hourly probabilities of 78 different crash types using generalized additive models. Using tensor product bases, we model non-linear relationships and combined effects of different meteorological parameters. We evaluate the increase in relative risk of different crash types in case of precipitation, sun glare and high wind speeds. The largest effect of snow is found in case of single-truck crashes, while rain has a larger effect on single-car crashes. Sun glare increases the probability of multi-car crashes, in particular at higher speed limits and in case of rear-end crashes. High wind speeds increase the probability of single-truck crashes and, for all vehicle types, the risk of crashes with objects blown on the road. A comparison of the predictive power of models with and without meteorological variables shows an improvement of scores of up to 24%, which makes the models suitable for applications in real-time traffic management or impact-based warning systems. These could be used by authorities to issue weather-dependent driving restrictions or situation-specific on-board warnings to improve road safety

Modeling hourly weather-related road traffic variations for different vehicle types in Germany

Weather has a substantial influence on people’s travel behavior. In this study we analyze if meteorological variables can improve predictions of hourly traffic counts at 1400 stations on federal roads and highways in Germany. Motorbikes, cars, vans and trucks are distinguished. It is evaluated in how far the mean squared error of Poisson regression models for hourly traffic counts is reduced by using precipitation, temperature, cloud cover and wind speed data. It is shown that in particular motorbike counts are strongly weather-dependent. On federal roads the mean squared error is reduced by up to 60% in models with meteorological predictor variables, when compared to models without meteorological variables. A detailed analysis of the models for motorbike counts reveals non-linear relationships between the meteorological variables and motorbike counts. Car counts are shown to be specifically sensitive to weather in touristic regions like seaside resorts and nature parks. The findings allow for several potential applications like improvements of route planning in navigation systems, implementations in traffic management systems, day-ahead planning of visitor numbers in touristic areas or the usage in road crash modelling

A classification algorithm for selective dynamical downscaling of precipitation extremes

High-resolution climate data O(1km) at the catchment scale can be of great value to both hydrological modellers and end users, in particular for the study of extreme precipitation. While dynamical downscaling with convection-permitting models is a valuable approach for producing quality high-resolution O(1km) data, its added value can often not be realized due to the prohibitive computational expense. Here we present a novel and flexible classification algorithm for discriminating between days with an elevated potential for extreme precipitation over a catchment and days without, so that dynamical downscaling to convection-permitting resolution can be selectively performed on high-risk days only, drastically reducing total computational expense compared to continuous simulations; the classification method can be applied to climate model data or reanalyses. Using observed precipitation and the corresponding synoptic-scale circulation patterns from reanalysis, characteristic extremal circulation patterns are identified for the catchment via a clustering algorithm. These extremal patterns serve as references against which days can be classified as potentially extreme, subject to additional tests of relevant meteorological predictors in the vicinity of the catchment. Applying the classification algorithm to reanalysis, the set of potential extreme days (PEDs) contains well below 10% of all days, though it includes essentially all extreme days; applying the algorithm to reanalysis-driven regional climate simulations over Europe (12km resolution) shows similar performance, and the subsequently dynamically downscaled simulations (2km resolution) well reproduce the observed precipitation statistics of the PEDs from the training period. Additional tests on continuous 12km resolution historical and future (RCP8.5) climate simulations, downscaled in 2km resolution time slices, show the algorithm again reducing the number of days to simulate by over 90% and performing consistently across climate regimes. The downscaling framework we propose represents a computationally inexpensive means of producing high-resolution climate data, focused on extreme precipitation, at the catchment scale, while still retaining the advantages of convection-permitting dynamical downscaling

Cell tracking of convective rainfall: sensitivity of climate-change signal to tracking algorithm and cell definition (Cell-TAO v1.0)

Lagrangian analysis of convective precipitation involves identifying convective cells (“objects”) and tracking them through space and time. The Lagrangian approach helps to gain insight into the physical properties and impacts of convective cells and, in particular, how these may respond to climate change. Lagrangian analysis requires both a fixed definition of what constitutes a convective object and a reliable tracking algorithm. Whether the climate-change signals of various object properties are sensitive to the choice of tracking algorithm or to how a convective object is defined has received little attention. Here we perform ensemble pseudo-global-warming experiments at a convection-permitting resolution to test this question. Using two conceptually different tracking algorithms, Lagrangian analysis is systematically repeated with different thresholds for defining a convective object, namely minimum values for object area, intensity and lifetime. It is found that the threshold criteria for identifying a convective object can have a strong and statistically significant impact on the magnitude of the climate-change signal, for all analysed object properties. The tracking method, meanwhile, has no impact on the climate-change signal as long as the precipitation data have a sufficiently high temporal resolution: in general, the lower the minimum permitted object size is, the higher the precipitation data's temporal resolution must be. For the case considered in our study, these insights reveal that irrespective of the tracking method, projected changes in the characteristics of convective rainfall vary considerably between cells of differing intensity, area and lifetime

Subhourly rainfall in a convection-permitting model

Convection-permitting models (CPMs)—the newest generation of high-resolution climate models—have been shown to greatly improve the representation of subdaily and hourly precipitation, in particular for extreme rainfall. Intense precipitation events, however, often occur on subhourly timescales. The distribution of subhourly precipitation, extreme or otherwise, during a rain event can furthermore have important knock-on effects on hydrological processes. Little is known about how well CPMs represent precipitation at the subhourly timescale, compared to the hourly. Here we perform multi-decadal CPM simulations centred over Catalonia and, comparing with a high temporal-resolution gauge network, find that the CPM simulates subhourly precipitation at least as well as hourly precipitation is simulated. While the CPM inherits a dry bias found in its parent model, across a range of diagnostics and aggregation times (5, 15, 30 and 60 min) we find no consistent evidence that the CPM precipitation bias worsens with shortening temporal aggregation. We furthermore show that the CPM excels in its representation of subhourly extremes, extending previous findings at the hourly timescale. Our findings support the use of CPMs for modelling subhourly rainfall and add confidence to CPM-based climate projections of future changes in subhourly precipitation, particularly for extremes